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O ! Abstract 

O 

In the Schwarzchild black hole spacetime, we show that chaotic motion can be triggered by 
the spin of a particle. Taking the spin of the particle as a perturbation and using the Melnikov 

hj ; 

\f\ . method, we find that the perturbed stable and unstable orbits are entangled with each other and 

<n : 

that illustrates the onset of chaotic behavior in the motion of the particle. 
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In astrophysics, the angular momentum or spin of a particle is one of the most important 
elements for consideration because nearly all astrophysical objects possess it. It may also be 
the crucial ingredient in the study of chaotic behavior as shown, for example, by Suzuki and 
Maeda [l| that the motions of spinning particles in the Schwarzschild black hole spacetime 
could be chaotic by analyzing the Lyapunov exponents |2j. Moreover, this is an important 
result due to its relevancy to the detection of gravitation waves from black hole coalescences 
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as it may jeopardize the use of matching templates in interpretating the upcoming 
data from the LIGO experiment. 

In this Letter we would like to use a different method, called the Melnikov technique 
[5J], to investigate the chaotic motions of spinning particles in the Schwarzschild black hole 
spacetime. In p| , we have successfully applied this technique to detect the onset of chaotic 
motions in the simple dynamical system of the one-dimensional Duffing hamiltonian with a 
spinning particle. In this method one consider homoclinic orbits emanating and terminating 
at the same unstable fixed point in an integrable system. When time-dependent perturba- 
tions are added to the system, the stable and the unstable orbits in the perturbed system 
would split in general. The Melnikov function, which is an integral evaluated along the un- 
perturbed homoclinic orbit, calculates the transversal distance between the perturbed stable 
and unstable orbits on the Poincare section. The isolated zeros in the Melnikov function 
indicated complicated entanglements between the two perturbed orbits and therefore the 
presence of chaotic behaviors. 

Here we first illustrate that there is a homoclinic orbit in the two-dimensional phase space 
of the radial motion of a spinless particle in the Schwarzschild spacetime. Suppose that the 
mass of the Schwarzchild black hole is M, and the line element is described in the usual 
Schwarzschild coordinates by pj 
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ds 2 = -fdt 2 + f- l dr 2 + r 2 (d6 2 + sin 

m = i - ^. (i) 

For a relativistic particle of mass m moving in a curved background described by the space- 
time metric g^, the action can be taken to be 

S[x] = -J g^Vdr. (2) 

x M (r) denotes the spacetime coordinates of the particle where r is the proper time, and we 



must in addition impose the mass shell constraint 

Pn = rng^x u , (3) 

where p^ is the momentum canonically conjugate to x M . Because g^ is stationary and 
axisysmmetric, the particle has two additional conserved quantities, the momenta conjugate 
to t and 4>, 

tp f dt 

E = -Pt = m J-f^ 

L = p^ = mr 2 sin 2 0—. (4) 

Let us restrict the particle motion, with a fixed angular momentum L, to lie in the 9 = 7r/2 
plane. Then the only effective degree of freedom is r, and we will find that there is exactly 
one homoclinic orbit in the reduced phase space (r,p r ) with appropriately chosen E and L. 
Substitute the expressions for E and L into Eq. ()3*|). a first-order equation of motion for 
r{r) can be obtained 

m2 © 2 + /(r) ( m2 + S) = B2 

\dr J r r 2 r 6 

This equation can be interpretated as describing the one-dimensional problem of a non- 
relativistic particle with energy E 2 — m 2 moving in the potential 

2Mm 2 L 2 2ML 2 
U -ff = — + ^ ~ -73-' ( 6 ) 

which has an attractive 1/r term and a repulsive 1/r 2 term as in the corresponding New- 
tonian system, plus an additional attractive 1/r 3 term which is a purely general relativistic 
effect. It is the latter which is responsible for the existence of the homoclinic orbit. 

It is convenient to change the dynamical variable to u = IMjr. The equation of motion 
in Eq. (jSJ) then becomes J8| 

,du^ AM 2 m 2 9 , AM 2 (E 2 -m 2 ) 
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FIG. 1: U eff (r), p = 0A. 

It is easy to see that C/ e // has extrema at real values of u, provided that 12M 2 m 2 / L 2 < 1, 
and they are located at 

1 



u s t 



;i-/9), 



U U n = -(!+/?)■ 



(8) 



where (3 = a/1 — l2M 2 m 2 / L 2 is the parameter in the problem to characterize the existence 
and the properties of the homoclinic orbit. u s t is the stable fixed point, while u un is the 
unstable one. Figure 1 shows the potential with f3 = 0.4. 

From Eq. (jBJ), we see that at the asymptotic value r = oo, the potential U e ff (oo) = . 
To obtain a homoclinic orbit, we must therefore require the condition U e ff(u un ) < 0, 



U eff {u un ) = -^(1 + /5) 2 (1 - 2(3) < 
0</3<|. 



(9) 



There could be a homoclinic orbit in this range of (3, and we will use the (3 = 0.4 to illustrate 
the chaotic behavior be in the homclinc orbit. Using the equations (3) and (5) we have the 
express for the momentum 



Vr 



mdr 1 

7^ = 7 



E 2 - f [m 2 + 
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(10) 



In Figure 2, we obtain a homoclinic orbit in the (r,p r ) reduced phase space, the parameters 
are m = 10" 6 , M = 1 and L = 3.78 x 10~ 6 for (3 = 0.4. 
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FIG. 2: p r (r), with m = 10" 6 ,M = 1,L = 3.78 x 1(T 6 . 

Next, we add spin to the particle on this homoclinic orbit to show that chaos may occur 
in the motion of a spinning particle around a black hole. The equations of motion of a 
spinning test particle in a general spacetime were first derived by Papapetrou £g . Those are 
a set of equations: 

— = u", 

(IT 

Dp" 1 

Dt ~ 2 "^ ' 

= p^-pW, (11) 

D/ Dt denotes a covariant derivative with respect to r. Here t> M , p", and S 1 ^ are the 4- 
velocity of the particle, the momentum, and the spin tensor, respectively. The multipole 
moments of the particle higher than the mass monopole and the spin dipole are ignored. It 
is called the pole-dipole approximation. 

There is a supplementary condition which gives a relation between v" and p" because p" 
is no longer parallel to v" in this case. We adopt the following condition [11], 

PrS*" = 0. (12) 

This condition is related to how to choose the center of mass in an extended body, and we 



can write down the relation between v^ and p^ explicitly, 
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(13) 



The equation shows that p^ is not parallel to v v in general due to the presence of the spin 
tensor. From Eq. (11) we get the equation of motion of the momentum p^, 

dp». — „ 1 



dr 



Vpu-^R^vs *, 



(14) 



which will be used in the analysis of the behavior of chaotic motions. 

Because the spacetime is static and spherically symmetric, there are two Killing vector 
fields, £& and (^ . The constants of motion are as given in [l( 
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C = C P1-1 ~ -C[i;vJ , 



E 



-a 



c*) 



-pt 



M qtr 
2 ' 



(15) 



Jz = C(j> = P<t> + r 2 sin 9 cos 



r sin 



E and J z are interpreted as the energy of particle and the z component of the total angular 
momentum, respectively. The x and y components of the total angular momentum are also 
conserved, 



J, 



•J n 



-pe sin <fi — p<f, cot 9 cos (j> 
—r sin <pS r + r 2 sin 2 9 cos 

•0 cos 4> — Pe), cot 6* sin <f> 
+r cos ^S" - + r 2 sin 2 # sin 



+ r sin # cos 9 cos 05 1 



+ r sin 9 cos sin 



tj>r 



(16) 



Because the background is spherically symmetric, we can choose the z axis to be in the 
direction of total angular momentum, then we have the constraint equations of the spin 
tensor, 

J 

— cot 9, 

_Pe_ 
r 



S re 



-\- J + ^h 

r V sin 9 



(17) 



From the constraint equation (12), we have the time-space components of the spin tensor, 



r = -U*+ "' 



rp t \ sin 9 

S te = — ( p r pe + -ps cot 9 
rp t \ r 

S» = ~— (jPr ~^ + -P0 COt e] • 

rp t \ sin 9 r J 



(18) 



Here we treat the spin tensor as a perturbation. Suppose that without these perturbative 
terms, the particle is originally in a homoclinic orbit on the equatorial plane 9 = ir/2. When 
spin terms are turned on, we write 

9 = --56. (19) 

Using 59 to qualify the degree of perturbation, we suppose that 

pe = 5p e + 0(59 2 ) 

P4> = J + 5 P(f) + 0(59 3 ) (20) 

where 5pg ~ 0(59) and 5p<f, ~ 0(59 2 ). With these we can expand the components of the 
spin tensor in Eqs. (17) and (18). 

S <t> = l59 + 0(59 2 ), 

S r8 = JM + 0(S9 2 ), 

r 
S* = -(J59 2 + 5p <p )+0(59 3 ), 

S tr = -^ (5p 2 e + J 2 59 2 + J5 P4> ) + 0(59 3 ), 

S te = ~ (p r Sp0 + —Se) +0(59 2 ), 
rh \ r J 

S^ = 4; (-JprSO 2 + -S Pe 59 - p r 5 P(P ) + 0(59 3 ). (21) 

rk, \ r J 

We have shown explicitly only the lowest order term in each spin tensor component in these 
equations. 

Similarly, we also expand the 4-velocity v M and the momentum of the particle p^ in 
Eqs. (13) and (14). In particular, we have 

f = —* + —i i 1 - 3J ^) m * + —M + °^- < 22 > 

dT mr 1 mr z \ m z r 6 J mr z 



Using this equation we change the parameter of the equations of motion from r to 0. The 
main equations of motion of the reduced phase space of the radial motion can then be written 

as 
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There are still the other three equations of motion from Eqs. (13) and (14), 



d56 

d(t> 



-jSp e + o(se 3 ), 
: ni ^^)se + o(se 2 ), 



-2M + 



3MJ 2 
m?r 2 



2M 



Pr 
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3M 



J 2 



m 2 r 2 



(25) 



5 Pe 56 



-^fpr^ + o^e 3 ). 

These equations are used to calculate the three perturbations 59, 5pg, and 5p^. We need 
these quantities as inputs in the main equations of motion, Eqs. (23) and (24). Since we want 
only the lowest order terms of these perturbations, we can take r and p r in these equations 
to be in their zeroth order, that is, corresponding to the unperturbed homoclinic orbit. 
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FIG. 3: 6$(<f>) 




FIG. 4: p e (cf>). 



Moreover, if we define e ~ 0(59), or 59 = e59, we can write 5po = e5pg and Sp^ = e 2x " 
with 59, 5pe, and Sp^ all of 0(1). Then the solutions with the initial conditions, 



59(0) = 1, 5p e (0) = Sp^O) = 0, 



(26) 



are plotted in Figs. 3-5, and the corresponding spin tensor components are shown in Figs. 6- 
8. We find them all to be oscillatory. Note that the choice of the initial conditions is 
quite arbitrary. We have checked the behaviors of those quantities for other sets of initial 
conditions like 59(0) = 0; 5p e (0) = 1; %(0) = 0, and 59(0) = 0; 5p e (0) = 0; %(0) = 1 . 
The results are also oscillatory as time evolves. 

We have illustrated that external oscillatory perturbative spin terms may induce chaotic 
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behaviors in simple systems [6j by calculating the Melnikov function p, |12[ . Now in the case 
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FIG. 7: 5 re (</>). 
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of a spinning particle in the Schwarzschild black hole spacetime, the Melnikov function, 

ME 2 7 - 1 ' ' x ~" 
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Here r and p r are again in their zeroth order corresponding to the unperturbed homoclinic 
orbit, with <fi — > + 0o, where 0o is used to parametrize the location along the unperturbed 
homoclinic orbit. Moreover, g\ and #2 involve only the oscillatory perturbative terms. The 
Melnikov function measures the transveral distance on the Poincare section between the 
perturbed stable and unstable orbits emanating from the unstable fixed point. The infinite 
number of discrete zeros of the Melnikov function as shown in Fig. 9 indicates that the two 
orbits entangle with each other, and this indicates the occurrence of chaotic behavior for 
the perturbed system. 

In addition, we have also calculated the Melnikov functions for the parameter j3 = 0.1, 
0.2, and 0.3. Similar results as shown in Fig. 9 with j3 = 0.4 are obtained. These values of j3 
satisfy the condition, < (3 < 0.5, for which unperturbed homoclinic orbits exist. We have 
thus proved the onset of chaotic motion for spinning particles near the homoclinic orbit in 
the Schwarzschild black hole spacetime , even with small spin terms. This is complementary 
to the results obtained by Suzuki and Maeda |l(, in which chaotic behavior is detected 
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FIG. 9: M(<f) ). 

numerically by calculating the Lyapunov exponent in the case of large spin terms. It is 
interesting to see if our method can be extended to consider other parts of the phase space 
and cases with finite spin terms. We plan to explore that in the future. 
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